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ABSTRACT 

We analyzed the results of non-linear resistive magnetohydrodynamical 
(MHD) simulations of jet formation to study the acceleration mechanism of ax- 
isymmetric, resistive MHD jets. The initial state is a constant angular momen- 
tum, polytropic torus threaded by weak uniform vertical magnetic fields. The 
time evolution of the torus is simulated by applying the CIP-MOCCT scheme ex- 
tended for resistive MHD equations. We carried out simulations up to 50 rotation 
period at the innermost radius of the disk created by accretion from the torus. 
The acceleration forces and the characteristics of resistive jets were studied by 
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computing forces acting on Lagrangian test particles. Since the angle between 
the rotation axis of the disk and magnetic field lines is smaller in resistive mod- 
els than in ideal MHD models, magnetocentrifugal acceleration is smaller. The 
effective potential along a magnetic field line has maximum around z ~ 0.5r in 
resistive models, where r$ is the radius where the density of the initial torus is 
maximum. Jets are launched after the disk material is lifted to this height by 
pressure gradient force. Even in this case, the main acceleration force around 
the slow magnetosonic point is the magnetocentrifugal force. The power of the 
resistive MHD jet is comparable to the mechanical energy liberated in the disk 
by mass accretion. Joule heating is not essential for the formation of jets. 

Subject headings: galaxies: jets — ISM: jets and outflows — accretion, accretion 
disks — diffusion — methods: numerical — MHD 



1. INTRODUCTION 

Magnetically driven mass outflows from accretion disks have been studied extensively. 
Blandford & Payne (1982) showed that a magneto-centrifugally driven cold outflow emanates 
from an accretion disk when the angle 9 between poloidal magnetic field lines threading the 
disk and the rotation axis of the disk is larger than 30°. Nonlinear magnetohydrodynamical 
(MHD) simulations of jet formation including accretion disks were first carried out by Uchida 
& Shibata (1985) and Shibata & Uchida (1986). Since the magnetically driven jet extracts 
angular momentum from the disk, the jet formation process enhances the accretion of the 
disk material. The back reaction of the jet formation on disk accretion and its relation to 
the magnetorotational instability (MRI ; Balbus & Hawley 1991) were discussed by Stone & 
Norman (1994), Matsumoto et al. (1996), and Kudoh, Matsumoto, & Shibata (2002). 

A key question which is often raised to the time dependent simulations of jet formation 
is whether the system approaches a steady state. Axisymmetric ideal MHD simulations 
of jet formation including an accretion disk show episodic outflows instead of approaching 
a quasi-steady state (e.g., Kuwabara et al. 2000). One successful approach to get steady 
jet by non-steady simulations is to treat the disk as time-independent boundary condition 
(e.g., Ustyugova et al. 1995; Romanova et al. 1997; Ouyed & Pudritz 1997; see Fendt & 
Cemeljic 2002, for resistive MHD simulations). Another approach is to include an accretion 
disk inside the computational box and assume magnetic diffusivity (Kuwabara et al. 2000; 
Casse & Keppens 2002, 2004). 

Kuwabara et al. (2000) carried out 2.5-dimensional axisymmetric resistive MHD simu- 
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lations starting from a rotating torus initially threaded by weak uniform vertical magnetic 
fields. They showed that the jet property drastically changes depending on the magnetic 
diffusivity, which they assumed to be uniform. When the magnetic diffusivity is small, mass 
accretion and jet formation takes place intermittently. On the other hand, in mildly diffusive 
disks, they showed by simulations for time scale about 50 rotation period at the innermost 
radius of the disk that both jets and accretion disks approach a quasi-steady state. Casse 
& Keppens (2002, 2004) extended this study to the case of geometrically thin disk initially 
threaded by equipartition (f3 = P gas /Pmag ~ 1) poloidal magnetic fields. By carrying out 
simulations assuming magnetic diffusivity localized inside the disk, they achieved near sta- 
tionary state. 

In order to study the acceleration mechanisms and energy transport of resistive MHD 
jets more quantitatively, we re-computed the simulation models adopted by Kuwabara et al. 
(2000) by applying the CIP-MOCCT scheme (Kudoh, Matsumoto, & Shibata 1998) which 
is more robust and accurate than the scheme we adopted in Kuwabara et al. (2000). 

In section 2, we describe the assumptions and numerical methods. Numerical results 
are presented in section 3. Discussions and conclusions are given in section 4. 



We solve two-dimensional nonlinear, time-dependent, compressible resistive MHD equa- 
tions in a cylindrical coordinate system (r, z) under the assumption of axisymmetry to in- 
vestigate the jet ejection from the accretion disk. The z-direction is parallel to the rotational 
axis of the accretion disk (see Fig. 1). The basic equations are 
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MODELS 



2.1. Assumptions and Basic Equations 



-f + V ■ Vp = -pV ■ V, 



(1) 



dV 



+ y-vy = — vp 

p 



dt 




(V x B) x B - Vip, 



(2) 



dB 



V x (V x B 



77V x B) , 



(3) 



dt 



— — h V ■ VP = -7PV 



V 



-4- 



-( 7 -l)£-(VxB) 2 , (4) 



where ip is the gravitational potential 



_ GM 

V- ( r 2 + z 2)l/2' W 

G is the gravitational constant, M is the mass of the central object, and rj is the resistivity 
which is assumed to be uniform. The other symbols have their usual meanings. In equation 
(4), we set r\ = rj or rj' = 0. The latter corresponds to the case when we neglect Joule 
heating (or equivalently, assume cooling which balances with the Joule heating). The units 
of length, velocity, time, and density are ro, Vko, to/Vko, and po, where tq is the radius where 
the density of the torus is maximum, Vko, and po are Keplerian rotation speed, and the 
density at (r, z) = (r , 0), respectively. In this normalization, we have two non-dimensional 
parameters: 

= (6) 
V 2 

FT — -M. (7\ 

^mg — T/2 , V I 

^KO 

where V s q = (7-P0/ Po) 1 ^ 2 , and Vao — [-^o/l^Po)] 1 ^ 2 are the sound speed, and Alfven speed at 
(ro, 0), respectively. Here, E th is the ratio of thermal energy to gravitational energy and E mg 
is the ratio of magnetic energy to gravitational energy. The normalized resistivity fj is defined 
as fj = rj/ (r VKo) and the magnetic Reynolds number at (r , 0) is defined as R m o = r V\o/v- 



2.2. Initial Condition 

We assume an equilibrium disk rotating around a central object surrounded by a hot 
corona (e.g., Matsumoto et al. 1996). The assumption of the existence of hot corona is a 
natural consequence of energy transfer from magnetically active disks. Such corona exists 
above galactic gas disks and solar photosphere (see e.g., Galeev et al. 1981, for accretion 
disk corona). Equilibrium solutions of a torus can be obtained under the following simpli- 
fying assumptions. Here, we adopted Newtonian analogue of Abramowicz's relativistic tori 
(Abramowicz, Jaroszynski, & Sikora 1978). The distributions of angular momentum is 

L = L r a . (8) 



We assume polytropic equation of state 

P = Kp 1+1/n . 



(9) 
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The density distribution of the torus is determined by 

GM 1 



2 „2a-2 



( r 2 + z 2)l/2 + 2 (l_ a ) L ° r 
P 

+(n + 1)— = constant. (10) 
P 

The mass distribution outside the torus is assumed to be that of the isothermal non-rotating 
high temperature halo surrounding the central object, 



P = Phexpja j- 2 + r ° 2)1/2 - 1 | 



(11) 



where a = iV^q/V^.. Here V sc and ph are the sound velocity and density in the halo at 
(0, r ), respectively. We assume that a = (L = constant), n = 3, 7 = 5/3, a = 1.0, 
Ph/ Po — 10~ 3 , E th = 5.0 x 10~ 2 , and E mg = 5.0 x 10~ 4 . The initial magnetic field is assumed 
to be uniform and parallel to the z-axis. 

The magnetic Reynolds number is defined by R m = XV a/ f] = R m o(^/ r o)(VA/VAo), where 
A = 2itVa/£1 is the characteristic scale of the magnetorotational instability, Va is the Alfven 
velocity, and Q is the angular velocity of the disk. It increases from inside the torus to the 
halo and it becomes R m ^> 1 in the halo, thus magnetic diffusion is not important there 
(Kuwabara et al. 2000). 



2.3. Numerical Methods and Boundary Conditions 

We solved the equations (l)-(4) by using the CIP-MOCCT method. The algorithm of 
the original CIP-MOCCT method and results of test simulations are described in Kudoh, 
Matsumoto, & Shibata (1998, 1999). In this method, the CIP scheme (Yabe & Aoki 1991) 
is used for hydrodynamical part and the MOCCT scheme (Stone & Norman 1992) is used to 
solve the induction equation (3) and to evaluate the Lorentz force terms. The basic equations 
(1), (2), and (4) are expressed in non-conservation form as 

^L + (V-V)f = S (12) 

where S is the source term. We solve the equation in two steps; advective step df /dt + 
(V ■ V)/ = and the source step df /dt = S. We revised the source step to ensure higher 
accuracy in space and time. In the MOCCT step, we modified the original scheme such that 
we include the resistive term in the induction equation and that characteristic equation of 
Alfven waves are solved by using the CIP scheme. 
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The size of simulation box is (r max x z max ) = (5.1r x 13.4r ), the number of grid 
points is (N r , N z ) = (200, 256), and the grid size is Ar = O.Olro when < r < 1 and 
Az = 0.01r when < z < 1. Otherwise the grid size increases with r and z. At r = 0, 
we assume p, P, V z , and B z are symmetric, while V r , V^, B r , and B^ are anti-symmetric. 
We computed only the upper half plane (z > 0) by assuming p, P, V r , V^, and B z are 
symmetric, while V z , B r , and B^ are anti-symmetric with respect to z — 0. The outer 
boundaries at r = r max and z = z max are free boundaries where waves can be transmitted. 
We softened the gravitational potential inside R — (r 2 + z 2 ) 1 / 2 = R m = 0.2 r to avoid 
the singularity at R — 0. As we show later, the disk material smoothly penetrates into 
the region R < Rm- R'm 

can be considered as the innermost radius of the disk. The unit 
time to = tq/Vkq ~ H-Rin/Vk,in where Vk,ui is the Keplerian rotation speed at r = R in , 
corresponds to 1.8 rotation period at r = i? in . 

We analyzed the results of three simulations (see Table 1). Model R is a mildly diffusive 
model (fj — fj' — 1.25 x 10~ 2 ). Model RC is the resistive model without Joule heating 
(fj = 1.25 x 10~ 2 , ff = 0). Model I is the non-diffusive model (rj = rj =0). Model R and 
Model / are the same as those reported by Kuwabara et al. (2000). 

To study the acceleration mechanism of the jet, we put Lagrangian test particles near 
the disk surface and computed the time evolutions of the location of these particles and 
evaluated the forces acting on each particle as schematically shown in Figure 1 (left). The 
initial positions of particles are selected such that they form main part of the jet. In model 
R, we put Lagrangian particles on a magnetic field line at r = 1.5 r , between z = 0.64 r and 
0.84 r . In model /, we put them on a magnetic field line at r = 0.8 r , between z = 0.31 tq 
and 0.46 r . The particle positions are updated by using fluid velocity V = (V r , V z ) as 
follows: r™ +1 = r™ + V r At, z™ +1 — z™ + V z At. Here, r p , and z p are the position of particles 
and n+ 1, and n show the time step whose interval is At. In Figure 1 (right), we schematically 
show the force f defined as the projection of poloidal force per mass f p to the direction of 
poloidal velocity vector V p of the particle; 

l Jpl I v pI 

When the sign of force /' is plus, the force accelerates the particle. On the other hand, 
when the sign of the force /' is minus, it decelerates the particle, respectively. By plotting 
the time variation of f , we can check the acceleration force along the streamline acting on 
particles. 
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3. 



NUMERICAL RESULTS 



3.1. 



Time Evolution of Jets 



Figure 2 (top panels) show the time evolution of temperature distribution (color scale), 
magnetic field lines (white curves) and velocity vectors (arrows) in model R. Numerical 
results are in good agreement with those reported by Kuwabara et al. (2000). As Kudoh 
et al. (1998) have already shown, numerical results obtained by the CIP-MOCCT scheme 
agree well with those obtained by the modified Lax-Wendroff scheme with artificial viscosity, 
except that numerical oscillations in low-/? regions are not prominent when the CIP-MOCCT 
method is used and that contact surfaces are sharply traced with the CIP-MOCCT method. 
Figure 2 (bottom panels) show the time evolution of model RC. In this model, resistivity 
is included in the induction equation but Joule heating term is not included in the energy 
equation. Numerical results indicate that Joule heating is not essential for the jet formation. 
The collimation of the jet is better in the model without Joule heating (model RC) because 
temperature of the jet is lower than that in the model with Joule heating (model R). 

Figure 3 (left) shows the density distribution at t — 25. High density ridge is formed 
near the outermost radius of the initial torus. This ridge is the contact surface between the 
disk material ejected from the disk and the ambient halo. The main part of the jet is inside 
this dense ridge. The density of the main part of the jet at z = 4ro is pj e t ~ 10~ 3 po but 
still larger than the halo density. Figure 3 (right) shows the volume rendered image of the 
density distribution (color scale) and the three-dimensional structure of magnetic field lines 
(solid curves). The magnetic field lines are highly twisted due to the rotation of the disk. 



Figure 4 shows the distribution of gas pressure (gray scale), magnetic field lines (white 
curves), and velocity vectors for model / (left), model R (middle), and model RC (right). 
In the ideal MHD model (Model /), gas pressure is small in the inner region (0.2 < r < 0.7 
and < z < 0.1), where magnetic pressure supports the disk. Mass accretion proceeds 
along the surface channel where the angular momentum of the infalling gas is magnetically 
extracted. Since the magnetic fields are frozen to the plasma, the mass accretion deforms the 
magnetic field lines. Magnetocentrifugal force accelerates the plasma along magnetic field 
lines which have sufficiently large angle from the rotation axis. 

In resistive models (model R and model RC), magnetic field lines are not deformed so 
much as those in the ideal MHD model because matter can traverse the magnetic field lines. 



3.2. 



Acceleration Force of Jets 
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Since the infalling matter loses less angular momentum than that in the ideal MHD model, 
centrifugally supported inner disk is formed. Due to adiabatic compression, gas pressure 
increases in the equatorial region. Gas pressure in the surface region of the disk is larger in 
model R than in model RC because of Joule heating. 

Figure 5 (right) shows the dependence of the inclination angle 8 of the poloidal magnetic 
field line in model R depicted by a white curve in the left panel. The angle from the 
rotation axis of the disk does not exceed 30°. This does not indicate that magnetocentrifugal 
acceleration is unimportant because fluid elements away from the equatorial plane can be 
accelerated by the magnetocentrifugal force even when < 30°. 

Figure 6 (left) shows the isocontours of temperature (gray scale contour), magnetic 
field lines (white curves), and velocity vectors (arrows) for model R at t — 13.1. Cold 
disk matter are accelerated along the magnetic field lines and form jets as we already 
showed in Kuwabara et al. (2000). Figure 6 (right) shows the distribution of pressure 
(gray scale contour), poloidal stream lines (white curves), and poloidal Lorentz force vectors 
(J x B) P /p = [(J x B) r , (J x B) z \ I p (arrows), where J = cV x B/(4n). The softened 
gravitational potential inside r = R in has small influence on the formation of jets emanating 
from the outer radius. The streamlines shown in Figure 6 (right) smoothly pass through 
R = i? in and slowly accretes to the central object. 

The poloidal Lorentz force is almost perpendicular to the poloidal stream lines in the 
launching region of the jet. This means that Lorentz force collimates the outflow toward 
the rotation axis but acceleration along the poloidal magnetic field lines is small. We would 
like to point out that inside the disk, the poloidal Lorentz force points toward the equatorial 
plane because the magnetic field mainly has +0 component and the electric current mainly 
has — r component. Thus the Lorentz force compresses the torus in the vertical direction. 
This force suppresses the outflow from this region. However, in the surface layer of the disk 
where the magnetic field lines change their direction from radial to vertical, the poloidal 
Lorentz force changes its direction and enables outflows. 

A white circle in Figure 6 (right) denotes the position of the test particle initially inside 
the torus at (r, z) = (1.5r , 0.73 r ) and later accelerated along the magnetic field line. For 
this test particle, poloidal Lorentz force almost has no contribution to the acceleration along 
the magnetic field line until t ~ 13.0. On the other hand, the test particle is in the region 
where pressure gradient force lifts the particle in the vertical direction. 

Figure 7 (left) shows the trajectories of the Lagrangian particles which are initially on 
a magnetic field line at r = 1.5 r in model R The initial position of particles is selected 
such that they form main part of the jet. Figure 7 (right) shows the time variation of the 
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forces per mass along the streamline exerted on the particle which is shown as open circles 
in Figure 7 (left). The forces accelerate the particle along its stream line when the sign is 
plus and decelerate the particle when the sign is minus. In Figure 7 (right), the curve F c 
shows the centrifugal force (v^/r)\\, the curve F g shows the gravity (— GM/r 2 )\\, the curve F p 
shows the pressure gradient force (—WP)\\/p, the curve Fl shows the poloidal Lorentz force 
per mass (J x B) p \\/p, and the curve F t shows the total force per mass obtained by summing 
the four forces. Here the subscript "||" means the component parallel or anti-parallel to the 
direction of the velocity vector of the particle. The two vertical broken lines show the time 
when V z changes from negative to positive (left line) and when V r changes from negative 
to positive (right line). Around t ~ 12.5, centrifugal force decelerates the radial inflow. 
The pressure gradient force contributes to turn the direction of motion toward the vertical 
direction. When the gravitational force becomes zero (t ~ 12.8), the pressure gradient 
force along the streamline becomes maximum, while the centrifugal force still decelerates 
the particle. At t — 13.16, centrifugal barrier finally turns the direction of the motion of 
the particle to +r direction. After the radial velocity becomes positive, centrifugal force and 
Lorentz force slowly accelerate the particle along the magnetic field line (t > 13.4). 

Figure 8 (left) shows the trajectories of particles which are initially located on a magnetic 
field line at r = 0.8 ro for model /. The initial position of particles is selected like that the 
particles form main part of the jet. Figure 8 (right) shows the forces along the streamline 
of the particle denoted by white circles in Figure 8 (left). Until t ~ 4.6, centrifugal force 
almost balances with the radial gravity (i.e., Keplerian rotation). Around t = 4.85 when the 
particle turns its direction to +r direction, the Lorentz force has the largest contribution to 
the acceleration of the particle. After t = 4.9, the summation of centrifugal force and the 
Lorentz force accelerates the particle. This result is consistent with that of Kudoh et al. 
(1998) and Kato, Kudoh, & Shibata (2002). Since the magnetic forces are much larger than 
those in model R (Fig. 7), acceleration is larger. 

Figure 9 shows the time dependence of r and z components of each force per mass 
acting on the Lagrangian particle shown as a non-filled circle in Figure 7 (left) for model R. 
The left panels show the r-component, and the right panels show the z-component. In the 
radial direction, centrifugal force exceeds the radial gravity around t ~ 12. Subsequently, 
radial component of the Lorentz force collimates the outflow (t ~ 14). In the vertical 
direction, the fluid element is lifted almost hydrostatically near the ejection time (t ~ 12.5) 
and subsequently accelerated by the Lorentz force (t ~ 14). It means that centrifugal barrier 
turns the direction of inflows. At the same time, vertical component of the pressure gradient 
force lifts the particle toward ^-direction and finally it is accelerated by the Lorentz force. 
Thus the radial inflow turns into the outflow along the magnetic field line. 
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Figure 10 (top) shows the magnetic pressure log 10 Pb and the gas pressure log 10 P g along 
a magnetic field line depicted in Figure 5 (left). The vertical dotted line in top, middle, and 
bottom shows the ejection point of the jet where V r — 0. At the ejection point, the gas 
pressure is dominant. On the other hand, the magnetic pressure becomes dominant in 
z > 0.35. Figure 10 (middle) shows the effective potential 

*- = -pri^ " \ aV ' (14) 

along a magnetic field line shown by a white curve in Figure 5 (left). The filled circle shows 
the slow magnetosonic point where V p = V s \ ow . Here V p is the poloidal speed and V^ ow is the 
slow magnetosonic speed defined as follows, 

1 



^slow 2 



V* + Vl-^(Vl + V?f-AVl A V s \ 



(15) 



where V p 2 A = B 2 p /(4np). 



The fluid elements in the region of d^ c s/dl < are accelerated by the magnetocentrifu- 
gal force and those in the region of d^ cS /dl > accrete in the case of ideal-MHD (Kudoh, 
Matsumoto, & Shibata 1998) where / is the line element of a magnetic field line. However, 
in the resistive model, the gas in the region d^/ e s/dl > is ejected because the gradient of 
gas pressure along a magnetic field line lifts the fluid elements and enable them to jump over 
the potential barrier. Subsequently, the gas is accelerated by magnetocentrifugal force. On 
the other hand, in ideal-MHD case, the ejection point locates at almost top of the effective 
potential (d^ e &/dl ~ 0) and the slow magnetosonic point almost coincides with this point 
(Fig. 10 bottom). 



3.3. Quasi Stationality of Outflow 

Figure 11 (top) shows the time evolution of fluxes of angular momentum and mass, 

Fl = 2 £ 5 2nr (rpV z V^ - "^^j dr, (16) 

/■2.5 

Fm = 2 2nrpV z dr, (17) 
Jo 

at z — 3r . The angular momentum and mass fluxes approach constant values around 
t = 25. Figure 11 (bottom) shows the distribution of the following quantities on a magnetic 
field line depicted in Figure 5. 

K = 4vrp(^, (18) 
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A = rV,- r -^, (19) 

n = *_^, (20) 

r Airpr 

S = log(^), (21) 

( 7 -l)p ^ 2 2 2 1 ' 

In ideal MHD case, these quantities should be constant along a magnetic field line (Ustyugova 
et al. 1999). The vertical dotted line shows the position where V r is zero (ejection point of the 
jet). The quantities are nearly constant even in the resistive model. These results indicate 
quasi-stationarity of the outflow. 



3.4. Powers of Jets 

To study the energetics of the jet formation and compare them with the results by Casse 
& Keppens (2004), we compute the energy liberated by accretion Pacc an d the power of jet 
Pjet defined as follows: 

-Pacc — Pmec + -Pent + Pvihd, (23) 

,0.3 (Y i x 

Pmec = -2 J 2nr pV r I— + *\ dz, (24) 

Pent = -2 I 2nr pV r ( ^—-) dz, (25) 

Jo \7-ipy 
r - 3 c 

Pmhd = -2J 2vrr — (E X B) r dz, (26) 

at r = 0.4 r , and 

Pjet = Pmec,j + Pent,j + Pvihd, j, (27) 
Pmec.j = 2 / 2nr pV z + A dr, (28) 



2.5 / r> 

7 P 







P ENT ,j = 2 / 2nrpV z [ ^37- ) dr, (29) 



f 2 - 5 c 

Pmhd,j = 2 J 2nr — (Ex B) z dr, (30) 
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at z = 3 r . 

Figure 12 shows the time evolution of the energy fluxes in the jet, a mechanical power 
Pmec.j, enthalpy flux Pent,j, Poynting flux Pmhd,j, and the ratio Ppovyj of each power to 
the total power Pjet at z — 3. The dominant energy flux at this height is the Poynting flux 
whose ratio to the total power is Pmhd.j ~ 70 %. The enthalpy flux and the mechanical flux 
contribute to Peng.j ~ 20 %, and Pmec.j ~ 10 %, respectively. 

Kudoh & Shibata (1997) showed that the dominant energy of a jet depends on the 
strength of magnetic field. When the poloidal component of magnetic field is B p oc r -2 , the 
fast magnetosonic point appears far from the Alfven point and the dominant energy of the 
jet is Poynting flux. In our simulations, initial magnetic field is uniform. In such models, 
the fast magnetosonic point locates far from the Alfven point (Kuwabara et al. 2000). 

Figure 13 (top) shows the magnetic flux P P E along a magnetic field line depicted in 
Figure 5, E (E oc r 2 ) is a cross section of the flux tube. Since P P E ~ constant, B p oc r~ 2 . 
Figure 13 (bottom) shows the total specific energy E = V p 2 /2 + + [7/(7 — 1)]P/ p + ip — 
rVL B^/ (Ati)\). Here V p /2 is the poloidal kinetic energy, and — rVL B^/ (AitX) is the Poynting 
flux divided by pV p where fi' = — V p B^jB p + V ( / } /r, and A' = pV p /B p . Figure 13 shows that 
Poynting flux is dominant. It is consistent with the case of B p oc r -2 in Kudoh & Shibata 
(1997). 

Figure 14 shows the time evolution of the mechanical power Pmec, enthalpy flux Pent, 
and Poynting flux Pmhd transported inward through r = 0.4 tq in the disk, and the ratio of 
the total power of the jet Pjet to the total power liberated by accretion Pacc- Numerical 
results of model R indicate that — Pent/Pviec ~ 0.4, Pmhd/Pmec ~ 0.1, and Pacc/Pmec ~ 
0.7. Meanwhile Pjet ~ Pa.cc- Thus, about 70% of the mechanical energy released by mass 
accretion powers the jet. 

4. SUMMARY AND DISCUSSION 

In this paper, we studied the acceleration mechanism of resistive MHD jets launched 
from accretion disks threaded by weak large scale poloidal magnetic fields. We re-computed 
the models we reported in Kuwabara et al. (2000) by applying the CIP-MOCCT scheme 
modified for resistive MHD equations. We carried out simulations for time scale about 50 
inner orbital time. Figure 15 shows the ejection mechanism of resistive MHD jets. In mildly 
resistive disks, the disk gas infalls without much losing its angular momentum because mag- 
netic field lines are less deformed than those in ideal MHD model. Thus, the matter accreting 
from the initial torus hits the centrifugal barrier and forms a high-pressure inner disk whose 
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pressure gradient force enables the accreting material to jump over the barrier of the effective 
potential. The material lifted up by pressure gradient force passes through the slow mag- 
netosonic point and is accelerated by the magnetocentrifugal force. Mass accretion/outflow 
takes place continuously and the system approaches a quasi-steady inflow-outflow state. 

The main acceleration force is the magnetocentrifugal force for both resistive and non- 
resistive models. The acceleration point by the magnetocentrifugal force depends on the 
resistivity and it is z ~ 0.6 r when fj ~ 0.01 and z ~ 0.25 r when fj = 0. The acceleration 
force for resistive model is only 25% of the acceleration for non-resistive case. We confirmed 
that Joule heating is not essential for the formation of jets. 

Through resistive MHD simulations treating the accretion disk as the fixed boundary, 
Fendt & Cemeljic (2002) showed that the jet velocity increases with increasing diffusivity. 
This result is consistent with the result of Kuwabara et al. (2000) that the mass outflow 
rate increases with the resistivity up to some critical value. The time scale approaching the 
quasi-steady state is shorter in our approach which includes the disk in simulation region 
because accretion of the disk material deforms the magnetic field lines. 

Casse & Keppens (2004) showed by MHD simulations assuming the magnetic diffusivity 
localized inside the disk that resistive MHD jets are formed from a thin disk threaded by 
global poloidal magnetic fields. Although the initial conditions and model parameters of our 
simulations are quite different from those of them, the density distribution and structures of 
magnetic field lines of the disk-jet system at the final quasi-steady state is similar (see Fig. 
2 of Casse & Keppens 2004). 

The energy flux of jets in our mildly resistive model is mainly transported by the Poynt- 
ing flux, while the mechanical flux dominates the outflow in Casse & Keppens (2004). Kudoh 
& Shibata (1997) showed that the dominant energy of a jet depends on the strength of mag- 
netic field. When the poloidal component of magnetic field is B p oc r~ 2 , the fast magnetosonic 
point appears far from the Alfven point and the dominant energy of the jet is Poynting flux. 
When B p oc r~( 2+a ) where a > 0, the fast magnetosonic point locates near the Alfven point 
and the dominant energy is the kinetic energy. In our model, B p is constant with radius 
at the initial state. Numerical results are consistent with those in Kudoh & Shibata (1997) 
when B p decreases slowly with radius. Since magnetic energy decreases faster in Casse & 
Keppens (2004), the fast magnetosonic point locates closer to the ejection point of the jet. 
Thus, in Casse & Keppens (2004), the jet is dominated by kinetic energy. 

The jet power (Pjet) is comparable to the energy released by mass accretion (Pacc)j 
which is about 70% of the released mechanical energy Pmec- As Casse & Keppens (2004) 
pointed out, in ADAF (Advection Dominated Accretion Flows) , Pmec ~ —Pent so that 
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-Pacc ~ 0. Thus, the resistive inflow-outflow configuration is different from ADAF. Our 
numerical results can be interpreted as magnetic analogue of advection dominated inflow- 
outflow solutions (ADIOS) which appears in hydrodynamical models with high phenomeno- 
logical viscosity (Blandford & Begelman 1999). 

We thank Dr. Seiichi Kato for helping us to extract forces acting on particles. This work 
was supported by the Japan Society for the Promotion of Science Japan-UK Cooperation 
Science Program (principal investigators: K. Shibata and N. O. Weiss), the Astronomical 
Data Analysis Center (ADAC) of the National Astronomical Observatory, Japan (NAOJ), 
ACT-JST of Japan Science and Technology Corporation, and National Science Council, 
Taiwan, Republic of China under the grants NSC-90-2811-M-008-020, NSC-91-2112-M-008- 
006. The numerical computations were carried out on VPP5000 at the ADAC of NAOJ 
under the projects rmhl7b, and rhn31b. 
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Fig. 1. — Left: Lagrangian test particles initially located near the surface of the accretion 
disk along a magnetic field line. Right: Schematic picture showing how to extract the force 
along the streamline. 
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Fig. 2. — Time evolution of temperature distribution (color contours), magnetic field lines 
(white curves), and velocity vectors (white arrows). Top panels: model R (with Joule heat- 
ing). Bottom panels: model RC (without Joule heating). 
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Fig. 3. — Left: Density distribution (color contours), velocity vectors (white arrows), and 
magnetic field lines (white curves) at t — 25. Yellow arrow shows the unit velocity. Right: 
Three-dimensional distribution of density (color) and magnetic field lines (blue curves). 



-19- 



Fig. 4. — Pressure distribution in Model / (left), Model R (center), and Model RC (right). 
White curves show the magnetic field lines, arrows show the velocity vectors and an arrow 
at upper right corner shows the unit velocity. 
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Fig. 5. — Left: Temperature distribution at t — 20.0 in model R (fj — 1.25 x 10~ 2 ). Right: 
Spatial variation of the angle 9 between the poloidal magnetic field line depicted in the left 
panel and the rotational axis as a function of z. Broken line shows the critical angle over 
which the magneto-centrifugal acceleration drives outflows. 
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Fig. 6. — Close up of the jet launching region of model R (fj — 1.25 x 1CT 2 ) at t = 13.1. 
Left: Gray scale contour shows the temperature distribution, the arrows show the poloidal 
velocity vectors (an arrow at upper right corner shows the reference velocity vector = I^ko) 
and white curves show poloidal magnetic field lines. White circle shows the position of the 
Lagrangian particle at this time. Right: Gray scale contour shows the pressure distribution, 
arrows show the poloidal component of the Lorentz force (an arrow at upper right corner 
shows the reference vector = 3GM/r 2 ) and white curves show poloidal stream lines. White 
circle shows the position of the Lagrangian particle at this time. 



Fig. 7. — Left: Trajectories of Lagrangian particles in model R (fj = 1.25 x 1CT 2 ). Dashed 
curves show the trajectories of Lagrangian test particles. Solid curve shows the initial disk 
surface. Circles are the location of test particles at times denoted in the figure. Right: Forces 
along a streamline of the particle shown by open circles in the left panel. Forces accelerate 
the particle when the sign is plus and decelerate the particle when the sign is minus. Radial 
velocity V r and vertical velocity V z change their sign at the time denoted by vertical broken 
lines. Curves F c , F g , F p , Fl, and F t show centrifugal force, gravity, pressure gradient force, 
poloidal Lorentz force, and total force. 
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Fig. 8. — Left: Trajectories of Lagrangian particles in model I (fj = 0.0). Dashed curves 
show the trajectories of test particles. Solid curve shows the initial disk surface. Circles are 
the location of test particles at times denoted in the figure. Right: Forces along a streamline 
of the particle shown by open circles in the left panel. Forces accelerate the particle when 
the sign is plus and decelerate the particle when the sign is minus. Radial velocity V r and 
vertical velocity V z change their sign at the time denoted by vertical broken lines. Curves F c , 
F g , F p , Fl, and F t show centrifugal force, gravity, pressure gradient force, poloidal Lorentz 
force, and total force. 
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Fig. 9. — Time evolution of r- and ^-components of each force acting on Lagrangian particle 
(non- filled circle) in Figure 7 (model R). 
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Fig. 10. — Top: Distribution of magnetic pressure Pb and gas pressure P g at t = 20 on 
a magnetic field line depicted in Figure 5. Middle: Distribution of Effective potential at 
t = 20 on a magnetic field line depicted in Figure 5. Bottom: Distribution of effective 
potential at t = 5.5 along a magnetic field line in main part of the ideal MHD jet (Model I). 
Vertical dotted line shows the ejection point of the jet (V r > 0). Filled circles show the slow 
magnetosonic point. 
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Fig. 11. — Top: Time evolution of fluxes of angular momentum Fl and mass Fm. Bottom: 
Distribution of variables, K = 4irp\ V p \/\B p \, A = rV^ — rB^/K, Q = V^/r — KB^/^npr), 
S = log(P/p 7 ), and E = [7/(7 - l)]P/p + ijj - fi 2 r 2 /2 + \/ 2 /2 + (V^/r - fi) 2 r 2 /2 at t = 20 
on a magnetic field line depicted in Figure 5. Vertical dotted line shows the ejection point 
of the jet (K > 0). 
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Fig. 12. — Time evolution of powers accelerating jet, mechanical power Pmec.J) enthalpy 
-Pent,j, Poynting flux Pmhd,j, and the ratio -Rpow.j of each power to total power. -Rmhd,j, 
-Rent, j, and -Rmec.j are the ratio of Poynting flux, enthalpy, and mechanical power to total 
power of jet. 
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Fig. 13. — Top: Magnetic flux B P T, at t — 20 along a magnetic field line depicted in Figure 
5. E(E oc r 2 ) is a cross section of the flux tube. Bottom: Specific energies at t = 20 along 
a magnetic field line in Figure 5. Solid curve shows the total energy E = V p 2 /2 + V^/2 + 
[7/(7 ~~ l)]-P/p + ^ ~ rQ' B^/(An\'), dashed curve shows the poloidal kinetic energy V p /2, 
and the dash-dotted curve shows the Poynting flux divided by pV p , —ril'B^/^irX'). 
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Fig. 14. — Time evolution of powers liberated by accretion, mechanical power Pmec> enthalpy 
-Pentj Poynting flux Pmhd transported inward through r = 0.4 r in the disk, and ratio of 
total power of jet to total power of accretion Pjet/|Pacc|; Pacc — -Pmec + -Pent + Pmhd- 
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Fig. 15. — Schematic picture of the ejection mechanism of the outflow in model R (fj — 
1.25 x 1CT 2 ). Gray region shows the high pressure inner torus created by the adiabatic 
compression of gases. 
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Table 1. Model Parameters 

Model fj RmO Outflow Type joule heating in energy eq. 

/ 0.0 oo Episodic no 

R 1.25 X 10~ 2 1.8 Quasi-steady yes 

RC... 1.25 X 10~ 2 1.8 Quasi-steady no 



Note. — In all models, E th = V s 2 /(-yV^ ) = 5 x 10" 2 , E mg = V* /V* Q = 
5 x 10~ 4 and p h /p = 10~ 3 
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